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ABSTRACT 

Gas disks of spiral galaxies can be described as clumpy accretion disks without a coupling 
of viscosity to the actual thermal state of t he gas. The model descr iption of a turbulent disk 



consisting of emerging and spreading clumps ([Vollmer fc Beckertll2003l ) contains free parameters, 



which can be constrained by observations of molecular gas, atomic gas and the star formation 
rate for individual galaxies. Radial profiles of 18 nearby spiral galaxies from THINGS, HER- 
ACLES, SINGS, and GALEX data are used to compare the observed star formation efficiency, 
molecular fraction, and velocity dispersion to the model. The observed radially decreasing ve- 
locity dispersion can be reproduced by the model. In the framework of this model the decrease 
in the inner disk is due to the stellar mass distribution which dominates the gravitational poten- 
tial. Introducing a radial break in the star formation efficiency into the model improves the fits 
significantly. This change in star formation regime is realized by replacing the free fall time in 
the prescription of the star formation rate with the molecule formation timescale. Depending on 
the star formation prescription, the break radius is located near the transition region between 
the molecular-gas-dominated and atomic-gas-dominated parts of the galactic disk or closer to the 
optical radius. It is found that only less massive galaxies (logM(M ) ^ 10)) can balance gas loss 
via star formation by radial gas accretion within the disk. These galaxies can thus access their 
gas reservoirs with large angular momentum. On the other hand, the star formation of massive 
galaxies is determined by the external gas mass accretion rate from a putative spherical halo of 
ionized gas or from satellite accretion. In the absence of this external accretion, star formation 
slowly exhausts the gas within the optical disk within the star formation timescale. 

Subject headings: galaxies: ISM — ISM: molecules 



1. Introduction (jBinnev et al. 2000), while the l ocal gas deple- 



tion time is w 2 Gyr (|Evandl2008l) . This suggests 



Gas accretion plays an important role for the that the gag con sumed by star formation must 

evolution of galactic disks. The star formation be rep i eil ished via external or internal accretion, 

rate in the solar neighborhood has been ap- Continuous addition of metal-poor gas would also 
proximately constant over the last ~ 10 Gyr 
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explain the discrepancy between the observed stel- 
lar metallicity distribution in the solar neighbor- 
hood and that pred icted by close d-box models of 
chemical evolution (|Tinslevlll98ll ). 

The need for accretion may also be seen from 
observations of other galaxies. Large spiral galax- 
ies like the Milky Way typically have star for- 
mation rates of a few M Q /yr and gas reservoirs 
of order 10 9 M Q in the star-forming part of the 
galaxy. Ov er these regions, ga s depletion times 
are 1-2 Gyr (jWong fc Blit j|2002f ) , so that for these 
galaxies to survive in their present configuration, 
roughly the same amount of mass consumed by 
star formation must arrive at the galaxy center 
via some type of accretion. This accretion can be 
external, infalling gas from the galactic halo, or 
internal, radial gas transport via spiral arms and 
gas viscosity. In the absence of external accretion 
the spiral galaxy might deplete its gas and become 
a lenticular galaxy. 

Turbulent galactic disks have been studied 
numerically and analy tically in recent years. 
IWada fc Normanl (|2007h used 3D, global hydro- 
dynamic simulations of the central part of a late 
type galaxy to study the density distribution of the 
inhomogeneous ISM. They found that the density 
distribution is well fitted by a single lognormal 



functi on over a wide density range. lElmegreenl 
(2002) first noticed that the Schmidt-Kennicutt 
law of the star formation rate can be reproduced 
if star formation occurs in dens e gas above a crit- 
i cal d ensity. The simulations of IWada fe Norman 
did not include feedback from supernova 
(SN) explosions. 



In a series of ar ticles , iKim fc Ostrikerl (l2007h. 



Shctt v fe Ostrikerl (|2008fl , and lKovama fe Ostriker 



( 20091 ) studied the 2D evolution of an inhomoge- 
neous ISM in a galactic disk without SN feedback. 



Kim fe Ostrikerl (|2007fl emphasize the impor- 



tance of gravity due to gas and stars for driving in- 
terstellar turbu lence. When strong feed back from 
SN is included (|Shettv fc Ostrikerll2008h . the frac- 
tion of dense gas and thus star formation is re- 
duced compared to the models of IKim fc Ostriker 
(|2007l ) and turbulence is driven by expanding 
shells in ove r dense regions. On the other hand, 
Agertz et al. I (l2009h . who studied gaseous galactic 
disks by means of 3D, high- resolution hydrody- 
namical simulations with and without SN feed- 
back, found that SN feedback is an important 



driver of turbulence in galaxies with star formation 
rates > 10~ 3 M Q yr _1 kpc~ 2 . This star formation 
rate is found to galactic radii ~ R25 in most of the 
nearby spiral galaxies investigated bv lLerov et al 



(120081). 

Shettv fc Ostrikerl |2008) note that the gas disk 



thickness is important for s etting the index of 



the S chmidt-Kennicutt law. iKovama fc Ostriker 
(2009) studied the relationships among pressure, 
vertical distribution of gas, and the fraction of 
dense gas in 2D hydrodynamic simulations of a 
starforming galactic disk with SN feedback. They 
found that vertical hydrostatic equilibrium gives 
a good estimate for the mean midplane pressure 
of the ISM. In their models they recover the ob- 
served ratio -Rmoi = Mh 2 /-Mhi if the gas surface 
density is proportional to the epicyclic frequency. 
iKoyama fc Ostrikerl (|2009i ) state that the empiri- 
cal result that -R mo i is proportional to the mean 
midplane pressure of the ISM implies that the 
epicyclic frequency, the gas surface density, and 
the star formation rate per unit volume are inter- 
dependent, most probably due to a galaxy evolu- 
tion with the Toomre parameter Q close to unity. 
If gravitational instability of a clumpy and tur- 
bulent gas disk is taken into account, turbulence 
drives the disc to a regime of transition between 
instability at small scales a nd stability in the cl as- 



sical sense (Toomre Q ~ l: lRomeo et all (|2010h ). 

Most recently, iKrumholz fc Burkertl (|2010h 
have developed an analytical model for the evo- 
lution of a thin disk of gas and stars with an 
arbitrary rotation curve that is kept in a state 
of marginal gravitational instability (Q ~ 1) and 
energy equilibrium due to the balance between 
energy released by accretion and energy lost due 
to decay of turbulence. Equili brium disks of this 
kind h a ve been investigated by Vollmer fc Beckertl 
(j2002h . IKrumholz fc Burkertl (|2010h showed that 
disks initially out of equilibrium evolve into it on 
timescales comparable to the orbital period if the 
external gas accretion rate is high. 

In this article w e use the analytical models of 



Vollmer fc Beckertl (|2003l hereafter VB03) to es- 



timate the radial mass accretion rate of a sam- 
ple of nearby spiral galaxies. The VB03 model 
treats galaxies as clumpy accretion disks using a 
simplified description of turbulence driven by su- 
pernova (SN) explosions. Because we focus on 
modeling the galactic disk within the optical ra- 
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dius, we neglect gravitational instabilities as en - 
ergy source for turbulence (see Agertz et al.l l2009) . 
A gas disk without star formation might generate 
turbulence via gravitational instabilities, but once 
stars form, we assume that SN dominate energy in- 
jection. The for malism of our equilibr i um m odel is 
close to that of iKrumholz fc Burkertl (|2010h . The 
model naturally links the physical properties of the 
galactic gas disk (surface density, molecular frac- 
tion, star formation rate) to the midplane pres- 
sure which mainly depends on the stellar surface 
density, epicyclic frequency, the Toomre parame- 
ter, and the gas mass ac cretion rate (similar to 
iKovama fc Ostriked 120091 ) . Furthermore, we as- 
sume the classical scale-independent Toomre sta- 
bility criterion. 

We compare these models to multiwavelength 
observa tions of 18 spiral ga laxies from t he THINGS 
sampl e (|Walter et al.112 008) compiled bv lLerov et al 



(120081) . From the comparison, we estimate the ra- 
dial mass accretion rate of each galaxy. 

The VB03 model considers the gas disk of a 
galaxy as a single turbulent medium. The dissipa- 
tion timescale of turbulent kinetic ener gy on galac- 
tic sc ales is about the crossing time ([Elmegreen 
2000), so that for a typical driving scale length 
of ~ 100 pc of the turbulent flow the dissipation 
time Td ~ 10 Myr . Therefore in order to maintain 
turbulence, an efficient, continuous, driving mech- 
anism is needed. Possible candidates for such a 
driving mechanism are: 



gravitati onal instabilities lik e spiral arms 
(see, e.g. iGomez fc Coxll2002l) or rotational 
s hear (see, e.g.. IWada fc Normanl 2001 ; 



Voll mer fc Beckert 2002; Kr umholz fc Burkert 
2010), 



Of these only SN explosions can balance the 
energy loss due to turbulent energy dissipa- 
tion if the driving lengt h scale is ~ 100 pc 
(|Mac Low fc Klessenll2004l) and the star formation 



rate is > 10" 3 Mgyr^kpc" 2 (|Agertz et alj |2009) 



In cases where the driving length scale is of the 
order of the disk thickness (~ 500 pc), the energy 
input due to gravitational instabilities can main- 
tain turbulence in the inte r stellar medium (ISM) 
(I Vollmer fc Beckerd l2002t IKrumholz fc Burkert 
2010). The VB03 model and this paper consider 
only turbulence driven by SN. 

In the VB03 model, SN-driven turbulence sets 
the disk structure and the disk structure deter- 
mines the star formation rate, which in turn de- 
termines the SN rate. As a result, the conver- 
sion of gas into stars is of prime importance to 
the model. Despite a large variety of observa- 
tions in different wavelengths, star formation in 
gal axies is not well understood (for recent reviews 



see 



Mac Low fc Klessenl 12004: Mc Kee fc Ostriker 



20071 ). ISM turbulence in turn probably con- 
trols the star formation proc ess to some de- 
gree (jMac Low fc Klessenl 120041) . but larg e scale 
gravitational inst abilities (iToomrel 1964 ), ther- 
mal instabilities ( Elmegreen fc Parravanol Il994t 



Wolfir e et al Jl995l) . and H2 formation (jKrumholz et al 



2009) may all be important. 



In the last decades it was thought that star 
formation is controlled by the interplay between 
self-gravitation and magnetic fields. However, the 
interstellar medium (ISM) is of turbulent nature. 
Its structur e is usually described as hierarchical 
(jScalol Il985h over length scales of several magni- 
tudes up to ^100 pc. The neutral phase of the 
ISM is not uniform but has a fractal structure 
(El megreen fc Falgarone 1996). 



magn eto-rotational instabilities (|Balbus fc Hawlevl |Vollmer fc Beckertl (|2003|) developed an analyt 



19911) . 



proto stellar outflows fe.g.. lWolf-Chase et al 



2000), 



stellar winds (e.g.. IVink et al 
UV radiation, 



20001) . 



ical model for clumpy accretion disks and included 
a simplified description of turbulence in the disk. 
In contrast to classical accretion disk theory (e.g., 
Pringlel[l98lh . VB03 do not use the "thermostat" 
mechanism, which implies a direct coupling be- 



tween the heat produced by viscous friction and 

the viscosity itself, which is assumed to be pro- 

I «ortional to the thermal sound speed. Thus, the 

st^er nova (SN) exptosio ns (see e.g.| Vollmer fc Beckerfl ^^ whkh ^ responsible for the gas heating5 



20031: lAgertz et al.l l20091 , 



external accretion ( Elmegre en fc Burkert 
20101 : iKlesse fc Hennebelldl2010l) . 



depends itself on the gas temperature. This leads 
to an equilibrium corresponding to a thermostat 
mechanism. Instead, VB03 use energy flux conser- 
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vation, where the energy flux provided by SNe is 
transferred through a turbulent cascade to smaller 
scales where it is dissipated. The SN energy flux 
is expected to be proportional to the local star 
formation rate. In particular, the star formation 
recipe takes the turbulent nature of the ISM into 
account. VB03 showed that it is possible to re- 
produce the the local and global gas properties of 
the Galaxy with this model. It was, however, not 
possible at that time to compare observed radial 
profile to the VB03 model. 

Here, we extend this co mparison to a sam ple of 



18 nearby spiral galaxies. iLerov et aD ( 20081 ) pre- 



sented radial profiles of the star formation rate and 
gas surface density (CO and Hi) of these galax- 
ies to study their star formation efficiencies. They 
showed that the SFE of H2 alone is nearly constant 
(H2 depletion time of 1.9 Gyr) at their 800 pc res- 
olution. Where the interstellar medium (ISM) is 
mostly Hi, however, the SFE decreases with in- 
creasing radius, a decline reasonably described by 
an exponential with scale length O.2R25-O.25R25. 

In this article we extend the formalism of VB03 
(Sec. [5]) by adding a break in the star formation 
rate (Sec. I2.1[) . Beyond this break the relevant 
timescale for star formation is no longer the free 
fall time of the most massive self-gravitating cloud 
at a given galactocentric radius, but the molecular 
formation timesc ale. In Sec. |Hand[5]the radial pro- 
files presented in lLerov et al.l (|2008l ) are compared 
and discussed to the extended VB03 model. We 
discuss the implications of our findings in Sec. [6] 
and give our conclusions in Sec. [7] 

2. The analytical model 

Since the model is described in detail in VB03, 
we summarize here only the basic idea and the 
resulting expressions for the relevant disk proper- 
ties. 

The model considers the warm, cold, and 
molecular phases of the ISM as a single, turbu- 
lent gas. We assume this gas to be in vertical 
hydrostatic equilibrium, with the midplane pres- 
sure balancing the weight of the gas and stellar 
disk. The gas is taken to be clumpy, so that the 
local density is enhanced relative to the average 
density of the disk. Using this local density, we 
calculate two timescales relevant to star forma- 
tion: the free fall of time of an individual clump 



and the characteristic timescale for H2 to form 
on grains. The longer of these (modulo an em- 
pirical scaling factor) is taken as the governing 
timescale for star formation. The star formation 
rate is used to calculate the rate of energy injec- 
tion by supernovae. This rate is related to the 
turbulent velocity dispersion (an observable) and 
the driving scale of turbulence. These quantities 
in turn provide estimates of the dumpiness of gas 
in the disk (i.e., the contrast between local and 
average density) and the rate at which viscosity 
moves matter inward. 

The model relies on several empirical calibra- 
tions: e.g., the relationship between star formation 
rate and energy injected into the ISM by super- 
novae, the H2 formation timescale (and its depen- 
dence on metallicity), and the turbulent dimension 
of the ISM (used to relate the driving length scale 
to the characteristic cloud size modulo a free pa- 
rameter) . As far as possible, these are drawn from 
observations of the Milky Way. 

We are left with three free parameters. First, 
there is an unknown scaling factor relating the 
driving length of turbulence to the size of gravita- 
tionally bound clumps, which we call 5. Second, 
the point at which the star formation timescale 
transitions from the free fall time to the H2 for- 
mation time is a priori unknown. Third, the mass 
accretion rate, M , which is related to the driving 
length and turbulent velocity, is a free parame- 
ter. In Section [¥] via comparison to the present 
sample of galaxies, we constrain 5. In practice, 
we also treat the Toomre Q parameter of the gas 
as a "semi-free" parameter, allowing it to change 
somewhat from the observed value. 

With these assumptions in hand, we may fit 
the model to a galaxy by comparing the observed 
kinematics, Hi, CO, and star formation rate pro- 
files to those predicted by the model. The result- 
ing fit yields the disk mass accretion rate and sev- 
eral other parameters that may be checked against 
expectations: the free fall time (or density) for 
the largest self-gravitating structures, the driving 
length of turbulence, and the approximate velocity 
of radial inflow. 

In the remainder of this section, we discuss our 
assumptions in slightly more detail, justify them 
via comparison to observation and theory, and 
note the physics that we neglect. 



4 



ISM: Following, e.g., Mac Low k Klessenl (|2004l ) 
we view the ISM as a single turbulent gas. In this 
picture, the warm, cold, and molecular phases 
of the ISM are a single entity. Locally, the ex- 
act phase of the gas may depend on the local 
pressure, metallicity, stellar radiation field, stel- 
lar winds, and shocks. Here we view these fac- 
tors as secondary, making a few simplifying as- 
sumptions. The equilibrium between the different 
phases of the ISM and the equilibrium between 
turbulence and star formation depends on three 
local timescales: the turbulent crossing time ^ urb , 
the molecule formation timescale t l mol , and the 
local free fall timescale t l s of a cloud. 

Supernova- driven Turbulence: First, we assume 
that the gas is turbulent, so that the turbulent 
velocity is the relevant one throughout the disk 
(making the exact temperature of the gas is largely 
irrelevant). We assume that this turbulence is 
driven by SNe and that they input their energy in 
turbulent eddies that have a characteristic length 
scale, /driv, and a characteristic velocity, t>t U rb- 
This driving length scale may be the character- 
istic length scale of a SN bubble, but it does not 
have to be so. It may be set by the interaction of 
multiple SN bubbles or of a SN with the surround- 
ing ISM. We note that based on simulations, the 
assumption of a single driving scale m ay be a sim- 
plification (|joung fc Mac Lowl Exit ) . The VB03 
model does not address the spatial inhomogene- 
ity of the turbulent driving nor the mechanics of 
turbulent driving and dissipation. It is assumed 
that the energy input rate into the ISM due to 
SNe, E$n, is cascaded to smaller scales without 
losses by turbulence. At scales smaller than the 
size of the largest selfgravitating clouds the energy 
is dissipated via cl oud contraction and star for ma- 
tion. We refer to Mac Low fe Klessenl ()2004h for 
a review of these topics. We limit our analytical 
model to the first energy sink which is the scale 
where the clouds become selfgravitating. 

We can connect the energy input into the ISM 
by SNe directly to the star formation rate. With 
the assumption of a constant initial mass function 
(IMF) independent of environment one can write 



Esn 



= £ S* — £ pJdiiv = Si/ 



c 'turb 

I 2 



(1) 



disk. The factor of proportionality £ relates the 
local SN energy input to the local star formation 
rate and is assumed to be independent of local 
conditions. £ is normalized using Galactic obser- 
vations by integrating over the Galactic disk and 
results in £ = 4.6 x 1CT 8 (pc/yr) 2 (see VB03). 
The adopted energy that is injected into the ISM 
is = 10 50 ergs bas ed on numerical studies 

bv lThornton et al.i (|l998l ). The final two parts of 
Equation [1] assume that stars form over a char- 
acteristic scale equal to the driving length and 
equate energy output from SNe with the energy 
transported by turbulence (see VB03). 

Star Formation in Molecular Clouds: Second, 
we assume that stars form out of gravitationally 
bound clouds. Where the timescale to form H2 is 
short, we take the local gravitational free fall time, 
given by 



iff = 



3tt 



32Gp c 



(2) 



to be the relevant timescale for star formation. 
Here G is the gravitational constant and p c \ the 
density of a single cloud. 

Cloud collapse, and thus star formation can 
only proceed if enough molecules form during the 
cloud collapse to allow the gas to continue cooling. 
Therefore, where the timescale for H2 formation is 
long (compared to the free fall time) , we view this 
as the relevant timescale for star formation. We 
take the characteristic time to form H2 out of H 
to be approximately 



(3) 



where AA is the unit surface element of the 



where a is a coeffici ent that depends on met allic- 
ity and temperature ( Draine fe Bertoldl [l996) and 
p c i is the density of a single cloud. Unlike some 
other recent numeri cal and theoretical treatments 
of H2 abundance bvlRobertson fc Kravtsov ( 20081 ) 
and Krum holz et al.l (l2009h . we take no account of 
the rate of destruction of H 2 via the UV radiation 
field. Instead, we assume H 2 to be destroyed via 
the collapse and subsequent star formation of a 
self-gravitating cloud. 

The coefficient of the molecular formation 
timescale a is assumed to be metallicity dependent 
(jTielens fc Hollenbachl Il985h . Because we admit 
external gas accretion, the metallicity of the star- 
forming ISM mainly depends on the ratio of ac- 
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cretion to star formation rate a. Small a < 1 lead 
to a metallicity derived from a closed box model, 
whereas in the case of a > 1 the metallicity equals 
the t rue stellar yield ytmo (jKoppen fc Edmunds! 
1999). For gas fractions higher than 0.04 the 
difference between the t wo solution s is les s than a 
factor of two. Moreover. iDalcanton ( 2007T ) showed 
that the effective yield y cS = Z gas /ln(l// gas ), 
where Z gas is the gas metallicity and / gas the gas 
fraction, for disk galaxies with a rotation velocity 
higher than 100 kms -1 is approximately constant, 
i.e., for these galaxies a closed box model can be 
applied. We thus feel confident to estimate the 
metallicity based on a closed box model using the 
gas fraction 0: 



a = 7.2xl0 7 x(log(- 



)) VMopc- 3 , (4) 



where £* is the stellar surface density and log 
refers to a natural logarithm. The calibration of 
the effective yield is based on the comparison of 
the observed O /H profiles with those of the closed 
box model (Sect. EJ. Adopting a stellar and gas 
surface density of E* = 40 M Q pc~ 2 and Ep- as = 
10 M (7)PC~ 2 at the solar radius of the Galaxy (jCox 
20Q0h yields a Q = 4.5 x 1 7 yrM^pc " 3 , which 



i s twi ce the value used by iHollenbach fc Tielens 
(Il997t i. 

This approach assumes that molecular clouds 
are relatively short lived, appearing and disap- 
pearing over roughly a free fall time (equivalently, 
by our construction, a turbulent crossing time); 
otherwise they might reach chemical equilibrium 
even when the H2 formation time is long com- 
pared to the free fall time. Accordingly, we es- 
timate the molecular fraction in the disk from 
the ratio of a cloud lifetime (i.e., the crossing 
or free-fall time) to the H2 formation time scale: 



Shi + Si-: 



= t 



turbAmol- 



fmo\ — 

Note that the molecular formation timescale, 
in particular, is fairly approximate — as noted 
above we have neglected radiation field and tem- 
perature dependences. Therefore we allow an ex- 
tra factor 7 when comparing the two to derive the 



iOnly three galaxies, IC 2574, NGC 4214, and NGC 2976, 
have rotation velocities smaller than 100 kms -1 . We cal- 
culated the models for these galaxies using up to 3 times 
smaller effective yields. The results where not significantly 
different compared to the models based on Eq. [4] with 
the exception that the molecular fraction decreased by the 
same factor as the effective yield. 



relevant timescale, with the transition between the 
two regimes at t l s = "ft l mol . We derive a value of 
7 = 0.4±0.3 for the VB03 star formation prescrip- 
tion and 7 = 0.12 ±0.06 for the star formati on pre- 
scription following iKrumholz fc Tanl (|2007j) . com- 
paring the model to the whole set of observations 
in Section [H 

Vertical Disk Structure: In the model, the disk 
scale height is determined unambiguously by the 
assumption of hy drostatic equilibr ium and the tur- 
bulent pressure (|Efmegreenlll989t ): 



Pturb = P^turb = 7T GS ( S + > ( 5 ) 

1 w disp 

where p is the average density, «t U rb the gas turbu- 
lent velocity in the disk, v*j isp the stellar vertical 
velocity dispersion, and £ the surface density of 
gas and stars. The stellar velocity dispersion is cal- 
culated by «j isp = y/2irGY l *H*, where the stellar 
vertical height is taken to be = Z y /7.3 with / 



being the stellar radial scale length (jKregel et al 



2002). We neglect thermal, cosmic ray, and mag- 
netic pressure. During the fitting procedure we 
realized that the local pressure equilibrium is of 
prime importance for the goodness of the fits. 

Treatment as an Accretion Disk: The turbulent 
motion of clouds is expected to redistribute an- 
gular momentum in the gas disk like an effective 
viscosity would do. This allows accretion of gas 
towards the center and makes it p ossible to trea t 
the disk as an accretion disk (e.g.. |Pringldll981l ). 



This gaseous turbulent accretion disk rotates in 
a given gravitational potential $ with an angu- 
lar velocity ft = \J R" 1 ^, where R is the disk 
radius. The disk has an effective turbulent viscos- 
ity that is responsible for mass accretion and out- 
ward angular momentum transport. In this case, 
the turbulent velocity is driven by SN explosions, 
which stir the disk and lead to viscous transport 
of angular momentum. In addition, star formation 
removes gas from the viscous evolution. Following 
Lin fc Pringld (|l987l ) , the evolution of the gas sur- 
face density is given by 

dY L _ 1 d { (d/dR)[isZR 3 (dn/dR)} \ ■ ■ 
~dt~~RdR\ {d/dR)(R 2 n) ) * + ° 

(6) 

where v is the gas disk viscosity, f2 the angular ve- 
locity, and S cxt is the external mass accretion rate. 
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In contrast to iLin fc Pringld (|l987l ) we assume a 
continuous and non-zero external gas mass accre- 
tion rate. By approximating d/dR ~ 1/-R, the 
global viscous evolution becomes 



as 

dt 



R 2 



(7) 



If the external and disk mass accretion rate keeps 
the combined Toomre parameter of the gas and 
stars smoothed over a few rotation periods to 
Qtot ~ 1, the gas surface density will only vary 
slowly with changes in the dark halo mass distri- 
bution (via fi) and the stellar disk structure, the 
gas loss due to star formation is balanced by exter- 
nal ac cretio n as suggested by iFraternali fc Binnev 
(|2008l ) and iMarinacci et all (|2010h . and the gas 
disk can be regarded as being stationary (dT,/dt = 
0). Indeed, local spiral galaxies s how a Qt.nt. 
not to o far away from unity (1.3-2.5: iLerov et al.1 
(2008)). For such a stationary gas disk, where star 
formation is balanced by external accretion, the 
local mass and momentum conservations yield: 



v M 
2?r 



(8) 



where M is the mass accretion rate within the 
disk. In the absence of external mass accretion, 
the gas disk can be assumed to be stationary as 
long as the star formation timescale i* exceeds the 
viscous timescale i„ = R 2 /v. For E cxt < £* and 
t* < t v Eq. |8]is not valid. In this case the gas disk 
is rapidly turned into stars within the gas con- 
sumption time (2 Gvr. lEvans (2008)). Since most 
spiral galaxies still have a significant amount of 
gas, we think that spiral galaxies are generally not 
in this state. Solving the time dependent Eq. [5] is 
beyond the scope of this work and we apply Eq. [5] 
Radial integration of Eq. [7] then gives 



dM 
~dt 



M — M* 4- M e - 



(9) 



where M is the disk gas mass. 

The viscosity is related to the driving length 
scale and characteristic velocity of the SN-driven 
turbulence by v = wturb^driv (VB03). Because 
the lifetime of a collapsing and starforming cloud 
(t l s < 4mb) * s smaller than the turnover time 
of the large-scale eddy (idrivAWb)) the turbulent 
and clumpy ISM can be treated as one entity for 
the viscosity description. 



dumpiness: A critical factor in the model is 
the relationship between the density of individual 
clouds, p c i, and the average density of the disk, 
p. It is the density of individual clouds that is 
relevant to the timescale for star formation. In this 
model, the two are related by the volume filling 
factor, $v so that p c i = 3>yV- 

Here p c \ refers to the density of the largest self- 
gravitating structures in the disk, so that for these 
structures the turbulent crossing time and gravi- 
tational free fall time are equal. The scale of such 
a cloud, Z c i is smaller than the driving length scale 
Idriv by a factor 5, which we do not know a priori. 

Shear due to differential galactic rotation could 
stabilize clouds, modifying the timescale for col- 
lapse. However, this effect is mainly important 
when the ratio of the cloud to disk surface den- 
sity is lower than the ratio of cloud to disk ve- 
locity dispersion, which is not the case over most 
of the disk in a typical spiral. Ty pical GMC sur- 
200 M pc" 2 (|Solomon et al 



face densities are 



1987T ). whereas disk surface densities only exceed 



100 M^pc 2 in th e very center of spiral galaxies 



([Lerov et al.l 12008). The ratio of velocity disper- 
sion being about 0.5, clouds at galactic radii larger 
than a few kpc are not stabilized by shear from 
galactic rotation. 

We can calculate the turbulent timescale for the 
cloud, t\ 



t[ ulh , for a fractal ISM: 



turb 



O 3 3 tdrivMurb 



(10) 



where D is the fractal dimension (see, e.g., 
Frischlll995l) of the ISM. We assume D = 2 for a 
compressible, self - gravitating fluid, which i s clos e 
to the findings of lElmegreen fc Falgarone) (|l996f ). 



Once S and thus t[ ulh are specified, we can solve 
for the density of the corresponding scale by set- 
ting t l s = t l tulh . The volume filling factor is then 
defined by comparing p c \ and p. 

Because the relationship between the driving 
length scale, Zdriv and l c \ is not known beforehand, 
we treat this as a free parameter. Once the vol- 
ume filling factor is known (from 5 or Z c i), we can 
calculate the local star formation rate, p\ 1 via 



P 



(11) 



sf 



where t l si is the local timescale for star forma- 
tion, either the free fall or molecular forma- 
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tion timescale, depending on local conditions. 
We assume n = 7 x 1 0~ 3 in consistency with 
Krum holz fc Tanl (|2007j ). Since in our model the 
lif etime of a cloud is the free-fall as su ggested 
bv lBallesteros-Paredes fc Hartmaim ( 2007 ) or the 
molecule formation timescale, this implies that 
during the cloud lifetime about 1 % of the cloud 
mass turns into stars. 

In VB03 a different star formation prescription 
is used: 



(12) 



The vertically integrated star formation rate in the 
inner disk where t\ f = t\ 



*turb — f turb IS 



p 



Iff 



urb 



(13) 



i.e. it is the mass flux density of the turbulent ISM 
into the regions of star formation. This alternative 
description yields fitting results that are as good 
as those based on Eq. [11] (in practice r\ ~ $y ~ 

10- 2 ). 

2.1. A Two Part Model 

VB03 assumed that t l s{ in EquationfTTlis always 
the local free fall time. As described above, here 
we consider two regimes, adopting instead the H2 
formation time when this becomes long compared 
to t l s . Therefore, the present model has two parts. 
We label the first, where iL = 
We label the second, where t l s 
disk. The change in regime should occur at the 
radius where t l s = 7t mol . Here 7 is a (free) scaling 
parameter that reflects a level of uncertainty in the 
application of the molecular formation timescale. 

An alternative interpretation of the change in 
star formation regime is that (i) the star forma- 
tion efficiency per free fall time or (ii) the vol- 
ume or area filling factor of star forming regions 
change from the inner to the outer disk . The latter 



t s the inner disk. 
: = t l mol the outer 



possib ility is motivated by the work of iBush et al 



( 20081 ) who show that gaseous spiral arms formed 
in the inner disk region can propagate into the 
outer gas disk creating overdensities where stars 
can form. Since these gaseous arms occupy a 
smaller fraction of the disk area with increasing 
galactocentric radius, the azimuthally averaged 
surface and area filling factors decrease. Having 



this picture in mind, one can write for the outer 
disk 



P* = V -T— = Vfmol^r = v4- (14) 



mol 



and 



p, = *v/- = <5>vf mo i4- = ^v-n 1 • (15) 

Vol ta tfl 



t\ t 1 



In this formulation the relevant local star forma- 
tion timescale is still the local free fall time of the 
clouds. 

2.2. Model calculations 

The VB03 model yields the following system of 
equations to describe a turbulent clumpy galactic 
accretion disk: 

^ ^turb^driv > 



M 



2tt ' 
Z = pH , 

Pturb = Kurb = o GE ( E + E 



Q 



2 

Vturb^ 

ttGE 



Vturb \ 

* ' 

disp 



y turb 



£* = $ v -^driv (VB03) or = V -^l dliv (KM05) , 



''SF 



*SF — 



3tt 



or 



32Gp c i Pel 

The meaning of the variables is given in Table [1] 
In the case of a selfgravitating gas disk (S ^> £*) 
or a dominating stellar disk (£* 3> S) the set 
of equations can be solved analytically (we use 
pv 2 twh = ttGY? and pv* nlb = ttGEE,). In Ap- 
pendix [A] we give the equations relating several 
observables («turb, S, SFE — £*/£, and / mo i and 
/mol = Sh 2 /(Sh 2 + Shi)) and other quantities 
of interest (Zdriv, v, <$>v) to our free parameters 
(M, S), assumed constants (£, 77), the approxi- 
mated metallicity (a), the "fundamental" observ- 
ables (stellar surface density £*, radius R and ro- 
tation f2) , and the "quasi-free" parameter Toomre 
Q of the gas (which we treat as free, but compare 
to observations). 
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For the global comparison between the observed 
and the model radial profiles we solve the set of 
equations given above numerically and we use 



L SF 



3tt 



32Gp c 



(16) 



3. Comparison with Observations 

We compare the model described above to ra- 
dial profiles of gas and star formation for 18 nearby 
spiral galaxies. 

3.1. Data 

The gal axies we compare t o are all part of the 
THINGS (IWalter et all l2008h survey. To carry 



out the comparison, we use radial profiles of the 
Hi velocity dispersio n measured via t he se cond 
moment presented by iTamburro et al. ( 2009 ) ; ro - 



tation curves derived by Ide Blok et al. ( 2008 1: 



radial profiles of atomic (Hi) and molecular gas 
(tr aced via CO) surfa ce density based on maps 
bv IWalter et all (120081) . iLerov et al.l (|2009l ). and 



iHelfer et al.l (j2003|); and profiles of star forma- 
tion r ate surface density derived by ILerov et"ail 
( 20081) from a comb ination of 24^m (SINGS, 



Kennicutt et~aH 120031) and FUV (the GALEX 
NGA. iGil de Paz et all [20071) intensity. We use 
stellar surface density e stimates d erived from 



SINGS 3.6^m imaging bv ILerov et al.l (|2008l ) as- 



suming a constant mass-to-light ratio. 

The profiles of atomic and molecular gas are 
combined into profiles of the molecular fraction, 
/ mo l = Eh 2 /(^h 2 + ^hi), and the profiles of star 
formation rate, Hi, and molecular gas are com- 
bined into profiles of star formation efficiency, 
SF E = Most of the profiles are drawn 

from lLerov et al. (2008) and details on methodol- 
ogy and data can be found there. 

The following uncertainties affect the data in 
addition to the uncertainty in the mean value over 
the azimuthal ring estimated from the RMS scat- 
ter in that ring: 

SFR and SFE: there is a ~ 50% uncertainty 
based on inter-comparison of different star forma- 
tion rate tracers. Half of this (~ 0.1 dex) is in- 
ternal to the galaxy, the other half (~ 0.1 dex) 
is galaxy-to-galaxy scatter (due to, e.g., the star 
formation history and dust properties). 



Hi: THINGS is expected to recover the true Hi 
flux with 10% accuracy. 

CO: The global CO-H2 conversion factor is un- 
certain by a factor of 2, which dwarfs any uncer- 
tainty in the calibration. An additional system- 
atic uncertainty is the dependence of the conver- 
sion factor on local conditions. The sample selec- 
tion somewhat minimizes this concern. The CO 
data have been masked before constructing pro- 
files, which can lead to artificially steep profiles at 
the edge of the CO emitting region. 

Velocity dispersion: the scatter in the disper- 
sion within a given ring is about 30%. Five of the 
18 galaxies are too inclined to derive reliable ve- 
locity dispersions. Several regions have very high 
velocity dispersions, > 20 km s _1 , which we expect 
are due at least partially to non-circular motions 
like bars, outflows, or streaming along spiral arms. 
We do not use these data in the comparison. 

3.2. Procedure 

The comparison is done in the following way: 

1. The Toomre parameter of the gas is calcu- 
lated using the observed Hi dispersion veloc- 
ity «disp, rotation velocity w ro t, and total gas 
surface density S gas : 



Q = 



^disp^rot 
7TGE gas i? 



(17) 



This approximate form of Q is used to avoid 
taking the (noisy) derivative of the observed 
rotation curve. In two cases, we adjust Q 
somewhat from the measured value: (i) in 
the inner part (R < 4 kpc) of the disk of 
NGC 3351 where we set Q = 8 and (ii) in 
the whole disk of NGC 2841 where we also 
assume Q = 8. For both galaxies this corre- 
sponds to an approximately constant veloc- 
ity dispersion of ~ 12-15 kms -1 . In cases 
without a measured velocity dispersion, we 
adopt a flat Q (NGC 2841, NGC 3351), an 
exponentially declining Q (NGC 3627), or 
a constant velocity dispersion of 15 kms -1 
(NGC 3521, NGC 7331). 

2. We solve the set of equations given in 
Sect. 12.21 using a grid for values of S, M, 
and 7. To account for the overall calibra- 
tion uncertainties (Sect. 13. ip . we calculate 
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the models on the <5-M-7-grid multiplying 
all points of Eh 2 by (0.63, 1., 1.58), tWb by 
(0.79,1., 1.26), and SFE by (0.71, 1., 1.41), 
allowing for all permutations. For each set 
of parameters the goodness g of the fit is 
determined for each observable profile: 



9 = 



1 N 

-Y 



obs 



V; 



model \ 2 



(18) 



where N is the number of points of a ra- 
dial profile, V^bs/modei are the observed and 
the model points of the profiles, and a is the 
uncert ainty of a point given by iLerov et al 



(12008b . At the end we calculate a total good- 
ness of the fit by taking the mean of the 
goodnesses of the individual profiles. This 
does not correspond to a traditional % 2 , be- 
cause the uncertainties of the points of a 
given profile are correlated. 

3. The final S and M are determined by the 
mean of all 5 and M yielding with min(g) < 
g < 1.1 mixi(g). We verified that these val- 
ues are not significantly different from those 
derived from the minimum of the total good- 
ness. In practice, the velocity dispersion pro- 
vides the strongest constraint on M, because 
it has a small error (it is shown in linear scale 
in Figs. [T]and[2]), followed by the molecular 
fraction and the SFE. Increasing M leads 
to an increasing velocity dispersion (Eq. IA1I 
and IA9|) and decreasing / mol (Eq. IA8I and 
IA161) and SFE (Eq. \M\ and lATSl Typi- 
cally, a variation in ¥ by a factor of 3 leads 
to changes in uturb, /mol> and SFE that are 
larger than the error bars (cf. Sect. [6]). 

4. Results 

We have applied the VB03 model using the 
KM05 (Eq. |rij and VB03 (Eq. |T2J star forma- 
ti qn descriptions to 1 8 spiral galaxies of the sample 
of lLerov etaH ( 20081) (Table [ 2 ]). Th e Hi rotation 
curves are from lde Blok et al.l (|2008 ). We realized 
that in the framework of the VB03 star formation 
prescription the free fall times of the largest self- 
gravitating clouds are sign ificantly larger than the 
collapse times derived by iTamburro et aL (2008) 
for NGC 0925, NGC 2403, and NGC 7793. Since 
the model fits are sensitive to the local vertical 



pressure equilibrium and the assumed correlation 
between the vertical and radial stellar scale height 
has a large scatter, we prefer a gravitational poten- 
tial dominated by the stellar disk for these galax- 
ies yielding loca l free f all times in agreement with 



Tamburro et al.1 (|2008l) . 



The resulting radial profiles are presented in 
Fig. [1] for the KM05 star formation prescription 
and in Fig. [2] for the VB03 star formation pre- 
scription. The order of the plots (from top left 
to bottom right) is: Hi rotation curve, observed 
and fitted stellar surface density profile E* (the 
assumed E* for the fit is shown as a solid line), 
model (dashed) and observed (solid) total gas den- 
sity profile E, Toomre Q parameter of stars+gas 
(dash-dotted) and of the gas derived from observa- 
tions (dotted) and assumed for the fit (solid), driv- 
ing scale length Zdriv, Hi velocity dispersion (solid) 
and model turbulent velocity dispersion (dashed) , 
observed (solid) and modeled (dashed) star for- 
mation efficiency, observed (solid) and modeled 
(dashed) molecular fraction, star formation (solid) 
and viscous (dashed) timescales, free fall timescale 
of the most massive selfgravitating gas clouds, ra- 
dial velocity of the gas within the disk, and the 
observed (solid) and model (dashed) metallicity. 
The observed metallicity profiles are d erived from 
12 + log(0/H) dMoustakas et al.ll2010h assuming a 
solar oxygen abundance of 12 + log(0/H) = 8.9. 
For the calculation of the combined To omre pa- 
rameter (stars+gas) Qtot 

(lRafikovll200l we lim- 
ited the scale of instabilities to 3 kpc. 

In agreement with Lerov et al. ( 20081 ) we find 
that Qtot is close to one for most of the galaxies. 
Only NGC 3521 and NGC 7331 show Q tot ~ 2 
and for NGC 2841 Q tot ~ 3. The Toomre param- 
eter of the gas ranges between 2 and 4 for most 
of the galaxies and these galaxies can thus be as- 
sumed in a quasi equilibrium state where the mass 
accretion (radially in the disk and external) bal- 
ances star formation to keep Qtot ~ 1- Moreover, 
for all galaxies with masses smaller than 10 10 M , 
the viscous timescale is smaller than or compa- 
rable to the star formation timescale. For these 
galaxies our assumption of a stationary accretion 
disk (Eq. [gj is justified without the need of a high 
external accretion rate. For more massive galax- 
ies Eq. [8] is a convenient parametrization which 
leads to a meaningful mean viscosity as discussed 
in Sect. [2] The parameters derived from the fitting 
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procedure are the disk mass accretion rate M, the 
scaling between the driving and dissipation length 
scale 5, and the break radius of the star formation 
efficiency -Rbrcak which are given in Table [3] for the 
KM05 and VB03 star formation prescriptions. 

We find acceptable fits for all galaxies (KM05 
and VB03), except for NGC 5194 (M 51). NGC 2976, 
NGC 4736, NGC 5194, NGC 3521, and NGC 5055 
show g > 1.5. For NGC 4736 we suspect a 
different mass accretion rate between the inner 
(R < 2 kpc) and outer disk. This would violate 
our assumption of stationarity, i.e. a radially con- 
stant mass accretion rate. For NGC 5055 (g ~ 4), 
the deviations between the model and the obser- 
vations occur in the central part of the galaxies 
where the velocity dispersion is overestimated. 
For NGC 2976 and NGC 3521 the disagreement 
between the model and observations occurs in 
the outer disk region. The KM05 star formation 
prescription provides better fits (smaller g) for 4 
galaxies (NGC 4214, NGC 4736, NGC 3351, and 
NGC 5194). The profiles of the other 14 galax- 
ies are better fitted by the VB03 star formation 
prescription. 

Whereas the KM05 star formation prescription 
yields a constant 6 — 1.9±0.5, 5 increases with 
galaxy mass in the VB03 star formation prescrip- 
tion. Introducing a 7 > leads to a significant 
improvement of the fit (about a two times smaller 
goodness). We find 7 = 0.10 ± 0.03 for the KM05 
star formation prescription. VB03 shows a bi- 
modal distribution of 7 with 5 galaxies having 
7 = 0.33 and 6 galaxies having 7 = 1.0. Over- 
all 7 = 0.65 ± 0.34 for the VB03 star forma- 
tion prescription. This translates into -Rbrcak/^* = 
3.0±1.1 and iZbrcak/^25 = 0.7±0.2 for the KM05 
star formation prescription and for massive galax- 
ies (log(M*) > 10) iJbreak/?* = 2.0 ± 1.0 and 
-Rbreak/^25 = 0.4±0.2 for the VB03 star formation 
prescription. The break radius of the VB03 star 
formation prescription is thus close to the tran- 
sition between a m ostly-Hi and a mostly-H2 ISM 
( Lerov et al. 2008), whereas it is closer to the op- 



tical radius in the KM05 star formation prescrip- 
tion. 

Total gas surface density: For most of the sam- 
ple galaxies the total gas surface density profile is 
well reproduced. In the framework of the KM05 
star formation prescription the total gas surface 
density is underestimated in the outer disk of 



IC 2574, NGC 7793, and NGC 0925; it is over- 
estimated in the outer disk of NGC 4736. In the 
framework of the VB03 star formation prescrip- 
tion the total gas surface density is somewhat un- 
derestimated in IC 2574, NGC 0925, and greatly 
underestimated in NGC 5194; it is overestimated 
in the outer disk of NGC 4736. All these differ- 
ences are within ~ 30%, except for NGC 5194 for 
the VB03 star formation prescription. 

Velocity dispersion: Velocity dispersion pro- 
files are available for 13 out of our 18 sample 
galaxies. We do not try to fit dispersion ve- 
locities higher than 20 kms -1 . These tend to 
appear in regions where geometry confuses the 
measurement (e.g., the central parts of inclined 
disks) or regions of very high SFR, where we 
worry that outflows or other non-disk structures 
render the measured velocity dispersion inappro- 
priate for comparison to the model. For exam- 
ples of o ur concerns, se e the p osition-velocity dia- 
grams of Ide Blok et al. (|2008l ) or compare the re- 
sults of gaussian fitting to moment methods treat- 
ing the same galaxy (e.g., iBoomsma et al. 20081 
Tamburro et al.l 120091 ). Both reveal the presence 



of asymmetric, multi-component Hi profiles in the 
central parts of some galaxies. We can reproduce 
the observed profiles within the error bars for all 
galaxies (KM05 and VB03 star formation prescrip- 
tions), except NGC 5194 with the VB03 star for- 
mation prescription. 

Star formation efficiency: The star formation 
efficiency in the inner disk of all galaxies is well re- 
produced by both models. The star formation effi- 
ciency of the outer disks of NGC 2976, NGC 0925, 
NGC 3198, NGC 5055, and NGC 7331 is better 
reproduced by the VB03 than by the KM05 star 
formation prescription. On the other hand, that of 
NGC 4214, NGC 6946, NGC 5194, and NGC 2841 
is better reproduced by the KM05 star formation 
prescription. NGC 2976 and NGC 4736 show very 
low surface density gas in this region leading to 
Q > 10. Since the star formation efficiency is pro- 
portional to Q" 1 in the outer disk, the mismatch 
between model and data might be due to our in- 
complete knowledge of Q in these regions. 

Molecular fraction: For 15 out of the 18 sam- 
ple galaxies CO observations are available. The 
VB03 star formation prescription reproduces the 
molecular fractions of almost all galaxies (except 
NGC 4736 and NGC 5194) roughly within the 
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error bars. There are deviations between the 
KM05 star formation prescription and observed 
molecular fractions for NGC 5194, NGC3521, 
NGC 5055. In NGC 3184, NGC 4736, NGC 6946, 
and NGC 5194 the decline of the observed molec- 
ular fraction at the outer edge of the CO distribu- 
tion is much steeper than predicted by the model. 
This is mostly due to sensitivity, because the CO 
data cubes were clipped at given S/N levels. 

5. Discussion 

Since the VB03 model is axisymmetric, it does 
not include non-axisymmetric structures like spi- 
ral arms and bars. Whereas spiral a rms do not 



chang e the given picture considerably ([Haan et al 



2009), the existence of bars do have a strong im- 
pact on radial gas flows. Despite the importance 
of a bar in the evolution of galactic disks, the fact 
that the VB03 model does not include their effect 
does not affect our conclusions. 

In the following we discuss aspects related to 
different radial profiles: 

Star formation law: Fig. |4] shows the star for- 
mation rate per unit area as a function of the 
total gas surface density (H1+H2) of our model. 
The distribution agrees well with observed distri- 
b utions from THIN GS and the literature (Fig. 15 
of lBigiel et al.ll2008l ). The absolute values and the 
shape, i.e. the 'knee', are reproduced by both the 
KM05 and VB03 models. 



Vel ocity dispersion: As shown in lTamburro et al 
(2009), the velocity dispersion of spiral galax- 
ies decreases with increasing galactocentric ra- 
dius. In the following, We note that the VB03 
model for a stellar disk dominated gravitational 
potential yields a radially declining velocity dis- 
persion, whereas a selfgravitating gas disk leads 
to a constant velocity dispersion in the inner disk 
(R < -Rbrcak) an d a declining velocity dispersion 
in the outer disk (R > i?brcak)- Thus, it is in prin- 
ciple possible to determine these different regimes 
from the radial behavior of the observed velocity 
dispersion. 

Molecular fraction: With the model radial pro- 
files we can investigate the dependence of the 
star formation rate per unit surface area on the 
molecular and total surface density (Fig. [3j). The 
correlation between the star formation rate and 
the molecular gas surface density is linear with a 



molecular depletion timescale of ~ 1.5 Gyr for the 
KM05 star formation prescription and ~ 2 Gyr 
for the VB03 star formation prescription. The 
correlation between the total gas surface den- 
sity and the star formation rate is steep where 
Hi is dominating (S < 20 M Q yr _1 ) and flat- 
tens for higher gas surface d ensities. The VB03 
mo del reproduc e s resu lts by iBigiel et~ (2008) 
and Lerov et al. ( 20081 ) who found (i) an average 
molecular gas depletion timescale of 2 Gyr and (ii) 
a critical gas surface density of 14 M pc -2 for the 
change between gas predominantly in molecular 
form and gas predominantly in atomic form. 

Driving length scale: The derived driving scale 
lengths are monotonically increasing with increas- 
ing radius. Typical values are 100-300 pc in 
the inner half of the optical disk. At the opti- 
cal radius the driving length scale is about 400- 
800 pc. The radial increase and the large val- 
ues at the optical radius are consiste nt with sizes 
of Hi shells observed in the G alaxy (lHeiledll984l 
McClure-Griffiths et alJliooih . 

Radial gas motions: Radial gas motions with 
velocity w rac i can be estimated with the mass ac- 
cretion rate M and t he gas surface density E gas 
(see e.g. IPringle 1981 ): 



M = 27riffig as (-i! ra d) 



(19) 



The radial profiles of the radial velocity are pre- 
sented in Figs. Q] and [5J Typically, we find ra- 
dial velocities smaller than 2 kms -1 . This is 



consis tent with the findings of iTrachternach et al 



(|2008l ) who found non-circular motions smaller 
than 5 kms -1 in these galaxies. Moreover, we ob- 
serve a general increase of radial velocities with de- 
creasing galactocentric radius. A strongly increas- 
ing Q in the outer disks of NGC 2976, NGC 4736, 
and NGC 3627 due to a strongly decreasing gas 
surface density leads to steeply increasing radial 
velocities (up to 10 kms -1 ) in these galaxies. This 
is again consistent with the profiles of n o n-circ ular 
motions derived by ITrachternach et al. ( 2008 ). 

Cloud free fall time estimates: The free-fall 
time only depends on the cloud density p c \ (Eq.[2]). 
Typical densities of giant molecular clouds are 
between 100-1000 cm -3 (jSolomon et all Il987t 
Hever et al. 20091) leading to f ree fa ll times of 



1.6-5.1 Myr. iTamburro et all (|2008f ) estimated 



a characteristic timescale for star formation in the 
spiral arms of disk galaxies, going from atomic 
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hydrogen (Hi) to dust-enshrouded massive stars. 
Their free fall time estimates vary between 1 and 
6 Gyr. The VB03 model yields the physical prop- 
erties of the most massive selfgravitating clouds 
at a given galactocentric radius via the local free- 
fall (ijj), turbulent (t[ uth ), and molecule formation 
(*moi) timescales. We recall that for selfgravitat- 
ing clouds t l s = 4urb- The radial profiles of the 
local free fall time for our sample galaxies are 
shown in Figs. [1] and [SJ They are calculated using 
Eq. [2] and IA101 For the comparison with the free 
fall timescales derived from observations we calcu- 
lated the mean and the standard deviation of the 
radial profiles within U < R < 21*. Our free fall 
timescales (Table [3]) are in good agreement with 
those derived from observations except for those 
of NGC 2403, NGC 0925, (both star formation 
prescriptions) and NGC 7793 (VB03 star forma- 
tion prescription) . By modifying Eq. |4] to better 
match the observed metallicity profile: 

E i y] 

a = 1.6x7.2xl0 7 x(log(— ^ — )) _1 yrM Q pc~ 3 , 

(20) 

we obtain the following parameters: NGC 7793: 
5 = 1, M=0.06 Moyr" 1 ; NGC 2403: 5 = 4, 
M=0.22 Moyr" 1 ; NGC 0925: (5 = 9, M=0.31 Moy^ 1 . 
With these parameters the free fall timescales of 
these three galaxies are in good agreement with 
expectations. We note that in this case 7 = 0.33 
for NGC 7793 and NGC 2403. 

6. How to sustain the star formation rate 

Can these galaxies sustain their star formation 
rates by radial transport of gas within the galactic 
disk? To answer this question one has to compare 
the local star formation rate and the local viscous 
timescale t v = R 2 jv using Eq. lA14l and lA22[ The 
local timescale comparison is presented in Figs. [1] 
and [5] The global comparison of the mean fraction 
< t v /t* > calculated over < R < R25 is shown 
in Fig. The KM05 and VB03 star formation 
prescriptions yield consistent results for this frac- 
tion within the galactic disks. 

The local viscous and star formation timescales 
are t v = R 2 jv = 2irY,R 2 / M and t., = £/£*. The 
mean fraction is then approximately < t u /t* >~< 
2ttY,*R 2 /M >~ M„/M. The important point for 
the global comparison is that the derived mass ac- 
cretion rates for almost all sample galaxies lie in 



the range between ~ 0.1 M Q yr _1 and 0.6 M^yr -1 , 
whereas the star formation rates vary between 
0.1 M^yr^ 1 and 3 Moyr -1 . In the following we 
show why the mass accretion rate shows this be- 
havior. Using v = fturb^driv in Eq. [T] leads to 



The expression for the mass accretion rate then 
becomes 

27rE 2 w 4 

M = 2tw£ = 27TUturbkrivE = r-^ . (22) 

££* 

Assuming typical values at R25/2, Uturb = 10 km s~ 
Zdriv = 100 pc, E = 10 M Q pc- 2 , yields a mass ac- 
cretion rate of M = 0.07 Moyr" 1 . For an error 
estimate we use the observational uncertainties 
given in Sec. 13.11 

Because of the large uncertainties associated 
with the molecular gas, we estimate the error 
of the mass accretion rate at galactic radii larger 
than -R25/2 where the gas is predominantly in 
atomic form. We assume E = 10 Mope -2 , 
AE = 1 M pc- 2 , E* = 5 x 10- 10 M pc- 2 yr~\ 
AE« = 2.5 x 10- 10 M Q pc- 2 yr~\ "tab = 
10 kms -1 , and Au tU rb = 3 kms" 1 . This leads 
to log^TXMoyr- 1 )) = -0.5 ± 0.5, i.e. we obtain 
an uncertainty of about a factor of 3. 

The inspection of Figs. [U [21 and[S]gives the con- 
sistent answer that only the less massive galaxies 
(log M* (Mg ) ;$ 10) can sustain the gas loss due 
to star formation by radial gas transport within 
the galactic disks. These galaxies can even have 
access to their gas reservoir beyond the optical ra- 
dius. On the other hand, the radial gas transport 
in the massive spiral galaxies might not be suffi- 
cient to balance the gas loss due to star forma- 
tion. This implies that, whereas a massive galaxy 
needs spherical infall from a putative gas halo or 
has to wait for infall with an angular momentum 
close to that of its disk to replenish its gas content, 
a less massive galaxy can live with large angular 
momentum accretion, because its mass accretion 
rate even at R > R25 is large enough for using 
this gas to sustain its star formation rate. The 
star formation rate of the massive galaxies is thus 
set by the amount of external accretion. In the 
absence of such an external gas accretion galax- 
ies will slowly consume their gas, the gas surface 
density will decrease, the Toomre parameter Q of 
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the gas will increase, and the star formation rate 
will decimal Examples in our sample might be 
NGC 3351 and NGC 2841. Given that most of 
the massive galaxies of our sample do not show 
this behavior suggests that these galaxies experi- 
ence mass accretion with rates comparable to their 
star formation rates (1-3 MQyr -1 ) from a putative 
spherical halo of ionized gas or from satellite ac- 
cretion leading to a temporarily enhanced mass 
accretion rate within the disk. 

7. Conclusions 

The theory of clumpy gas disks (VB03) provides 
analytic expressions for large-scale and small-scale 
properties of galactic gas disks. Large-scale prop- 
erties considered are the gas surface density, den- 
sity, disk height, turbulent driving length scale, ve- 
locity dispersion, gas viscosity, volume filling fac- 
tor, and molecular fraction. Small-scale proper- 
ties are the mass, size, density, turbulent, free-fall, 
molecular formation timescales of the most mas- 
sive selfgravitating gas clouds. These quantities 
depend on the stellar surface density, the angular 
velocity f2, the disk radius R, and 4 free parame- 
ters, which are the Toomre parameter Q of the gas, 
the mass accretion rate M, the ratio S between the 
driving length scale of turbulence and the cloud 
size, and the radius, at which the local star for- 
mation timescale is no longer the cloud free-fall 
timescale, but the molecule formation timescale. 
We determine these free parameters using three 
independent measurements of the radial profiles 
of the (i) neutral gas (Hi), molecular gas (CO), 
and star formation rate (FUV + 24 u rn). A sam- 
ple of 18 mostly spiral galaxies from iLerov et al.l 
(|2008l ) is used in the analysis. Based on the simul- 
taneous VB03 model fitting of the radial profiles of 
the total gas surface density, velocity dispersion, 
star formation efficiency, and molecular fraction, 
we conclude that 

1. the model star formation efficiency is very 
sensitive to the description of local pressure 
equilibrium in the disk midplane (Eq. [S]) ; 

2. the fits of all radial profiles are acceptable for 
all galaxies except NGC 5194 (M 51). The 



2 This decline is so mewhat slowe d down by gas replenishment 
from dying stars (Gislcr 1979). 



model-derived free-fall timescales of selfgrav- 
itating clouds are in good agreement with 
expectations from observations; 

3. the observed radially decreasing ga s veloc- 
ity dispersion (|Tamburro et al.l [2009) can be 
reproduced by the model. In the frame- 
work of the VB03 model the decrease in the 
inner disk is due to the stellar mass dis- 
tribution which dominates the gravitational 
potential. A selfgravitating gas disk yields 
a constant velocity dispersion in the inner 
disk, whereas it leads to a radially decreas- 
ing velocity dispersion in the outer disk. It 
might be thus possible to identify the differ- 
ent regimes from the radial behavior of the 
gas velocity dispersion; 

4. Introducing a change in star formation 
regime into the model improves the fits sig- 
nificantly. This change is realized by replac- 
ing the free fall time in the prescription of 
the star formation rate with the molecule 
formation timescale; 

5. depending on the star formation prescrip- 
tion, the best-fit break between regimes 
in the model is located near the tran- 
sition region between the molecular-gas- 
dominated and atomic-gas-dominated parts 
of the galactic disk or closer to the optical 
radius; 

6. the viscous timescale is smaller than or com- 
parable to the star formation timescale for 
galaxies less massive than 10 10 M , whereas 
it is much higher for the massive galaxies; 

7. as a consequence less massive galaxies can 
balance the gas loss due to star formation 
by radial gas inflow within the galactic disk. 
In this way these galaxies can even use the 
gas reservoir outside the optical radius. This 
is impossible for massive galaxies. The star 
formation rate of massive galaxies is deter- 
mined by the external gas mass accretion 
rate. Massive galaxies depend thus on exter- 
nal infall with an angular momentum close 
to that of the disk, whereas less massive 
galaxies can use large angular momentum 
gas located beyond the optical radius for star 
formation. 
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A. Model equations 

A.l. Star formation recipe according to Krumh olz &: Ta n (2007) (Eq.HH 



The following equations are appropriate for the inner disk regime where the local free fall time is the 
limiting timescale for star formation and E* 3> £■ For the actual numbers we assume M — 0.2 M©yr _1 , Q = 
2, E* = 100 Moyr" 1 , S = 2, t] = 7 x lO" 3 , a = 7.2 x 10 7 x (log(^S)) _1 yrM Q pc~ 3 = 2 x 10 7 yrM Q p C - 3 , 
7 = 0.2, and Q, = 1.8 x 10~ 8 yr" 1 . 

Vtmh = l.08G*M*Q*6*Z*ri*'£]!tl-i = 19 kms^ 1 , (Al) 

0v = 0.50G%Mh-%Q-ich~^**n* = 0.01 , (A2) 

£ = 0.50G-fM*Q"t5^*r?*s|o§ = 16 M pc~ 2 , (A3) 
S, = riS^Q- 1 ^ = 1.3 x 10~ 8 M yr _1 pc -2 , (A4) 

l dl . iv = 0.43G3MlQij-^-^-iE*^f2-t = 140 pc , (A5) 

v = 0A6G%M§Qi5-i£-ir)-iY,**n-i = 2.8 x 1CT 3 pcV" 1 , (A6) 

SFE = 1.99 GiM-*Q-*5*^*r?*E|fi* = 7.7 x lO" 10 yr" 1 , (A7) 

/moi = O^G^M-aQ-^^T^Ef fi^aT 1 = 0.80 . (A8) 

The following equations are appropriate for the inner disk regime where the local free fall time is the 
limiting timescale for star formation and E>E,: 





(A9) 


(j) V = 0.62G%Mh-¥£-i v -i = 0.02 , 


(A10) 


E = 0.28G~$M*Q- X 8$£*rjin = 7 M Q pc~ 2 , 


(All) 


£* = 0.28G-*A/*g- 2 (5^577§^ 2 = 9.2 x 10~ 10 Moyr^pc -2 , 


(A12) 


/driv = 0.66G^Mf QS-i^-^rj-iQ- 1 = 350 pc , 


(A13) 


v = O.blGHl^QS'iC^rj-in- 1 = 6.0 x 10~ 4 pcV" 1 , 


(A14) 


SFE = rjSQ^fl = 1.3 x 10~ 10 yr^ 1 , 


(A15) 


/ mol = 0.39 G~^ Q~ l 8^ rfi SlaT 1 = 0.21 . 


(A16) 



Note that the turbulent velocity and the volume filling factor are constant and thus independent of the 
galactocentric radius. 

The following equations are appropriate for a selfgravitating outer disk where the molecular formation 
timescale is the relevant timescale for star formation and E 3> E»: 

Vtuth = O.MG-iQ-i8£ir)i("/a)-ini = 12 kms -1 , (A17) 

0v = 10.6G 5 Af 2 Q 3 (5" 8 r 3 f?" 3 a 3 ^ _3 = 0.01 , (A18) 

E = 0.17G-5Q-I ^ijjia-lnl = 7 M pc" 2 , (A19) 

£* = 0.02 G- 5 M ^Q-^^^jfa^n 5 = 1.4 x 10" 9 Moyr^pc -2 , (A20) 

Z driv = 1.7 G 2 MG a 8~ 2 £~ 1 'n~ 1 aDr 2 = 330 pc , (A21) 

v = 0.92 G^MQiS-^-^-ia^Q-i = 3.0 x 10~ 3 pc 2 yr _1 , (A22) 

SFE = 0.09 G-i M~ l Q-i 8 h ^rfi (7a) - 5^ = 1.4 x 10~ 10 yr" 1 , (A23) 

/ mo i = QMG-^M^Q-iS^^^a-ini = 0.23 . (A24) 



15 



A. 2. Star formation according to VB03 (Eq. H2|) 

The following equations are appropriate for the inner disk regime where the local free fall time is the 



limiting timescale for star formation and E* ^> E. For the actual numbers we assume M = 0.2 M Q yr 1 , 
Q = 2, £» = 100 M yr-\ (5=2,a=7.2x 10 7 x (log(^S))~ 1 yrM p C - 3 = 2 x 10 7 yrM Q p C - 3 , 7 = 0.2, 
and ft = 1.8 x 10~ 8 yr" 1 . 

v turb = imG&M&Q&S-^^Epn-Ts = 19 kms" 1 , (A25) 

cjjy = 0.73GAM*^-iiQ-w^-T5S^*fi^ = 9 x 10" 3 , (A26) 

E = 0.33 G~ A_/vf Aq-t3(5-tS£t3 £F£lT3 = 12 M Q pc~ 2 , (A27) 

£* =0.73GAmAq-tI5-A^-toS?Qt5 = 1.5 x 10" 8 Meyr^pc^ 2 , (A28) 

Zdriv = 0.48G*M^Qw5ra^-HE;*fi-A = 140 pc , (A29) 

v = 0.49GwMT§QT5^^-wS7^fi-T5 = 2.8 x 10" 3 pc 2 yr" 1 , (A30) 

12 • 2 3 5 8 10 o — 

SFE = 2.23G^MT3Q-T3,5-T3^-T3S* 13 rjT3 = 1.3 x 10~ 9 yr" 1 , (A31) 

/ mo i = OmG-^M-^Q-^S^^Zpn-ha- 1 =0.88 . (A32) 

The following equations are appropriate for the inner disk regime where the local free fall time is the 
limiting timescale for star formation and E > E»: 

Wturb = 0.82 G*M*<r tt£tt = 13 kms" 1 , (A33) 

4.4 1R Q 

</>v = 0.81 GttM"<T Trf-TT = 0.01 , (A34) 

E = 0.26 G _ *M*Q -1 <rnf irft = 8 M pc~ 2 , (A35) 

£* = 0.21 G~* ji/TTQ-^-TT^-Trft 2 = 1.4 x 10~ 9 Moyr^pc- 2 , (A36) 

Z driv = 0.71 G*Af AQ^-*^ 1 = 290 pc , (A37) 

v = 0.61 G^M^QS^C^^ 1 = 4.1 x 10~ 3 pc 2 yr" 1 , (A38) 

SFE = 0.81 G^M^Q-^rAf-^n = 1.9 x 10~ 10 yr" 1 , (A39) 

/ mo i = 0.34G-TTA'/-AQ-i ( 5Tr^TT^ Q! -i = 0.27 . (A40) 



Note that the turbulent velocity and the volume filling factor are constant and thus independent of the 
galactocentric radius. 

The following equations are appropriate for a selfgravitating outer disk where the molecular formation 



timescale is the relevant timescale for star formation and E>E,: 

vturb = 0.71GBM3Q-tr£tr( 7a )-!ns = 13 kms^ 1 , (A41) 

0v = 1.95G3Af5Q!,r 2 £-! (7a) t ft" I = 9 x 10~ 3 , (A42) 

E = 0.23G-SM 4Q-i£t?( 7a )-i?ft§ = 8 M pc" 2 , (A43) 

E* =0.09G-iM5Q-T£-i( 7a )-!ftT = 1.7 x 10~ 9 Moyr'^c" 2 , (A44) 

Zdriv = 0.93 G*M*Qi ( 7 a)-ft _ - = 270 pc , (A45) 

v = 0.69G»MiQw£~»(7a)«ft~» = 3.9 x 10~ 3 pcV" 1 , (A46) 

SFE = 0.39G-iM^Q-T^-i( 7Q: )-iftT = 2.1 x 10~ 10 yr" 1 , (A47) 

/ mo] = 0.22G"^M"3Q-T^I 7 ~l a - J rr2T = 0.29 . (A48) 
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A. 3. Critical angular velocity 

Here we give the expressions for the critical angular velocity f2 C rit where tg = ^ft mol (S 3> £*)• For 
the actual numbers we assume M = 0.2 M yr _1 , Q = 2, £* = 100 Mgyr" 1 , S = 2, r\ = 7 x 10~ 3 , 
a = 7.2 x 10 7 x (log(S^))" 1 yrM Q pc- 3 = 2 x 10 7 y r Mgp c" 3 , 7 = 0.2, and O = 1.8 x 10~ 8 yr" 1 . 

Star formation recipe according to Krumholz fc Tan ( 2007f ) fEq. ITT|): 



n ciit = 2.6 G ■* Ms (J-s £s r)Sj a Q = 1.7 x 10" 8 yr" 1 . (A49) 
Star formation according to VB03 (Eg. [T2"|) : 

fi crit = 2.9GTTM *5-*^-A 7 aQ = 1.3 x 10~ 8 yr" 1 . (A50) 

Assuming a constant rotation velocity of v Iot = 200 kms -1 , we obtain critical radii of 12.4 kpc and 
16.2 kpc, respectively. 
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Fig. 1. — Comparison between modeled (KM05) and observed radial profiles. Form upper left to lower right: 
(i) Hi rotation curve, (ii) observed and fitted stellar surface density profile E* (the assumed E* for the fit 
is shown as a solid line), (iii) model (dashed) and observed (solid) total gas density profile E, (iv) Toomre 
Q parameter of stars+gas (dash-dotted) and of the gas derived from observations (dotted) and assumed for 
the fit (solid), (v) driving scale length IdHv, (vi) Hi velocity dispersion (solid) and model turbulent velocity 
dispersion (dashed), (vii) observed (solid) and modeled (dashed) star formation efficiency, (vii) observed 
(solid) and modeled (dashed) molecular fraction, (ix) star formation (solid) and viscous (dashed) timescales, 
(x) free fall timescale of the most massive selfgravitating gas clouds, (xi) radial velocity of the gas within the 
disk, and (xii) the obs erved (solid) and model (dashed) metallicity. The observed metallicity profiles is derived 



from 12 + log(0/H) (jMoustakas et al.ll2010h assumin g a solar oxygen abundance of 12 + log(Q/H) = 8.9. 



The t wo solid lines correspond to the calibrations of iKobulnickv fc Kewlevl (|2004D and iPilvugin fc Thuan 
(|2005j ) with an additional offset of 0.6 dex. 
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Fig. 2. — Comparison between modeled (VB03) and observed radial profiles. Form upper left to lower right: 
(i) Hi rotation curve, (ii) observed and fitted stellar surface density profile E* (the assumed E* for the fit 
is shown as a solid line), (iii) model (dashed) and observed (solid) total gas density profile E, (iv) Toomre 
Q parameter of stars+gas (dash-dotted) and of the gas derived from observations (dotted) and assumed for 
the fit (solid), (v) driving scale length IdHv, (vi) Hi velocity dispersion (solid) and model turbulent velocity 
dispersion (dashed), (vii) observed (solid) and modeled (dashed) star formation efficiency, (vii) observed 
(solid) and modeled (dashed) molecular fraction, (ix) star formation (solid) and viscous (dashed) timescales, 
(x) free fall timescale of the most massive selfgravitating gas clouds, (xi) radial velocity of the gas within the 
disk, and (xii) the obs erved (solid) and model (dashed) metallicity. The observed metallicity profiles is derived 



from 12 + log(0/H) (jMoustakas et al.ll2010h assumin g a solar oxygen abundance of 12 + log(Q/H) = 8.9. 



The t wo solid lines correspond to the calibrations of iKobulnickv fc Kewlevl (|2004D and iPilvugin fc Thuan 
(|2005j ) with an additional offset of 0.6 dex. 
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Fig. 3. — Correlation between the star formation rate per unit surface area and the molecular gas surface 
density. The solid line correspond to a molecular gas depletion rate of 2 Gyr. Left panel: KM05 star 
formation prescription. Right panel: VB03 star formation prescription. 
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Fig. 4. — Correlation between the star formation rate per unit surface area and the total gas surface density. 
The solid line correspond to a molecular gas depletion rate of 2 Gyr. Left panel: KM05 star formation 
prescription. Right panel: VB03 star formation prescription. 
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Fig. 5. — Mean fraction between the model-derived mass accretion timescale and the observed star formation 
timescale for Z* < R < -R25. The error bars are based on a factor of 3 uncertainty of the mass accretion rate 
(see Sect. [5]). Left panel: KM05 star formation prescription. Right panel: VB03 star formation prescription. 
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Tabic 1: Model Parameters. 



Parameter 


Unit 


Explanation 


G = 5 x icr 15 


pc 3 yr" 1 Mn i 


gravitation constant 


K 


yr -1 


epicyclic frequency 


Q 




Toomrc parameter 


R 


pc 


galactocentric radius 


H 


pc 


thickness of the gas disk 


-H* 


pc 


thickness of the stellar disk 


-^break 


pc 


break radius between star formation regimes 




pc 


cloud size 


Vrot 


pcyr -1 


rotation velocity 


= V ro t/R 


yr" 1 


angular velocity 






volume filling factor 


$A = *v H/l c i 




area filling factor 


p 


M pc- 3 


disk midplane gas density 


Pel = P/$V 


M pc- 3 


cloud density 


P* 


M Q pc- 3 yr- 1 


star formation rate 


s 


M pc- 2 


gas surface density 




M pc- 2 


stellar surface density 




M Q pc- 2 yr- 1 


star formation rate 


£ = 4.6 x 1CT 8 


p C 2 yr -2 


constant relating SN energy input to SF 


M 


Moyr" 1 


disk mass accretion rate 


«turb 


pcyr -1 


gas turbulent velocity dispersion 


«rad 


pcyr" 1 


gas radial velocity 


U disp 


pcyr" 1 


stellar vertical velocity dispersion 


V 


9 — 1 

pc yr 


viscosity 


/mol = Sh 2 /(ShI " 


I-Sh 2 ) 


molecular fraction 


a 


yrM Q pc 3 


constant of molecule formation timescale 


^driv 


pc 


turbulent driving length scale 


^diss 


pc 


turbulent dissipation length scale 


(5 




scaling between driving and dissipation length scale 




yr" 1 


star formation efficiency 




yr 


cloud free fall timescale at size I 


Hurb 


yr 


cloud turbulent timescale at size I 


''mol 


yr 


cloud molecule formation timescale at size I 






star formation efficiency per free fall time 


7 




ig f = 7i mol at the break radius 
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Table 2: Properties of Sample Galaxies (from Lcroy et al. 2008). 



Galaxy 


Dist. 




P.A. 


Morph. 


M B 


-R25 


log M* 


SFR 






(Mpc) 


n 


(°) 




(mag) 


(kpc) 


(M ) 


(Moyr- 1 ) 


(kpc) 


IC 2574 


4.0 


53 


56 


Irr 


-18.0 


7.5 


8.7 


0.07 


2.1 


NGC 4214 


2.9 


44 


65 


Irr 


-17.4 


2.9 


8.8 


0.11 


0.7 


NGC 2976 


3.6 


65 


335 


Sc 


-17.8 


3.8 


9.1 


0.09 


0.9 


NGC 7793 


3.9 


50 


290 


Scd 


-18.7 


6.0 


9.5 


0.24 


1.3 


NGC 2403 


3.2 


63 


124 


SBc 


-19.4 


7.3 


9.7 


0.38 


1.6 


NGC 0925 


9.2 


66 


287 


SBcd 


-20.0 


14.2 


9.9 


0.56 


4.1 


NGC 0628 


7.3 


7 


20 


Sc 


-20.0 


10.4 


10.1 


0.81 


2.3 


NGC 3198 


13.8 


72 


215 


SBc 


-20.7 


13.0 


10.1 


0.93 


3.2 


NGC 3184 


11.1 


16 


179 


SBc 


-19.9 


11.9 


10.3 


0.90 


2.4 


NGC 4736 


4.7 


41 


296 


Sab 


-20.0 


5.3 


10.3 


0.48 


1.1 


NGC 3351 


10.1 


41 


192 


SBb 


-19.7 


10.6 


10.4 


0.94 


2.2 


NGC 6946 


5.9 


33 


243 


SBc 


-20.9 


9.8 


10.5 


3.24 


2.5 


NGC 3627 


9.3 


62 


173 


SBb 


-20.8 


13.9 


10.6 


2.22 


2.8 


NGC 5194 


8.0 


20 


172 


SBc 


-21.1 


9.0 


10.6 


3.13 


2.8 


NGC 3521 


10.7 


73 


340 


SBbc 


-20.9 


12.9 


10.7 


2.10 


2.9 


NGC 2841 


14.1 


74 


153 


Sb 


-21.2 


14.2 


10.8 


0.74 


4.0 


NGC 5055 


10.1 


59 


102 


Sbc 


-20.6 


17.4 


10.8 


2.12 


3.2 


NGC 7331 


14.7 


76 


168 


SAb 


-21.7 


19.6 


10.9 


2.99 


3.3 
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Table 3: Model results for the KM05 and VB03 star formation prescriptions. 



Galaxy 


M 


8 


7 


-f^brcak 


CH 2 


c v diBp 


CSFE 


9 


< t s > 




(Meyr- 1 ) 






(kpc) 










(Myr) 


KM05 


IC2574 


0.22 


2.0 


0.07 


1.1 


0.63 


0.79 


0.71 


0.98 


29.6 


NGC4214 


0.08 


1.3 


0.11 


1.3 


0.63 


0.79 


0.71 


0.15 


9.5 


NGC2976 


0.31 


3.2 






0.63 


1.00 


1.00 


2.59 


4.1 


NGC7793 


0.06 


1.4 


0.11 


3.9 


0.63 


0.79 


0.71 


0.69 


6.0 


NGC2403 


0.22 


1.8 


0.07 


4.9 


1.00 


1.00 


0.71 


0.31 


8.8 


NGC0925 


0.11 


2.0 


0.07 


9.1 


0.63 


0.79 


0.71 


0.97 


12.5 


NGC0628 


0.08 


1.6 


0.13 


7.6 


0.63 


1.00 


1.00 


0.95 


4.7 


NGC3198 


0.11 


1.8 


0.07 


11.0 


0.63 


0.79 


1.00 


0.94 


7.2 


NGC3184 


0.11 


2.0 


0.10 


10.0 


0.63 


1.00 


1.41 


1.24 


3.8 


NGC4736 


0.89 


2.0 






0.63 


1.26 


1.41 


2.34 


7.7 


NGC3351 


0.06 


1.6 


0.10 


8.6 


0.63 


1.00 


0.71 


0.79 


5.8 


NGC6946 


0.32 


2.5 






0.63 


1.00 


1.00 


1.47 


3.0 


NGC3627 


0.22 


2.0 






0.63 


0.79 


1.00 


0.47 


4.4 


NGC5194 


1.26 


1.0 






0.63 


0.79 


0.71 


3.29 


22.1 


NGC3521 


0.22 


1.8 






0.63 


0.79 


1.41 


5.82 


4.9 


NGC2841 


0.08 


2.0 


0.17 


11.3 


0.63 


0.79 


0.71 


0.42 


6.0 


NGC5055 


0.32 


2.0 


0.07 


12.5 


0.63 


0.79 


1.41 


7.83 


5.9 


NGC7331 


0.22 


2.0 


0.09 


13.9 


0.63 


0.79 


1.41 


1.37 


4.4 


VB03 


IC 2574 


0.16 


8.9 


1.00 


0.00 


0.63 


0.79 


0.71 


0.21 


5.2 


NGC 4214 


0.06 


1.0 


1.00 


0.00 


0.63 


0.79 


0.71 


0.21 


11.2 


NGC 2976 


0.32 


3.5 






0.63 


1.00 


0.71 


2.44 


4.1 


NGC 7793 


0.16 


1.0 


1.00 


0.00 


0.63 


1.00 


0.71 


0.33 


12.5 


NGC 2403 


0.22 


1.3 


1.00 


0.00 


0.63 


1.00 


0.71 


0.11 


14.2 


NGC 0925 


0.16 


1.0 


0.33 


0.00 


0.63 


0.79 


0.71 


0.41 


32.0 


NGC 0628 


0.08 


1.8 


1.00 


0.00 


0.63 


1.00 


1.00 


1.01 


4.6 


NGC 3198 


0.11 


2.5 


0.33 


4.3 


1.00 


0.79 


1.00 


0.93 


5.1 


NGC 3184 


0.16 


4.0 


0.33 


6.7 


1.58 


1.00 


1.41 


1.31 


2.3 


NGC 4736 


0.63 


7.9 






1.58 


1.26 


1.41 


2.59 


1.7 


NGC 3351 


0.06 


1.8 


0.33 


2.7 


0.63 


1.00 


0.71 


0.84 


5.1 


NGC 6946 


0.45 


7.9 






1.58 


1.00 


1.00 


1.13 


1.2 


NGC 3627 


0.32 


8.9 






1.58 


0.79 


0.71 


0.40 


1.3 


NGC 5194 


0.06 


8.9 






1.58 


0.79 


0.71 


4.30 


0.8 


NGC 3521 


0.32 


8.9 


0.25 


11.7 


1.58 


0.79 


1.41 


4.22 


1.1 


NGC 2841 


0.11 


3.5 


0.33 


5.1 


1.00 


0.79 


0.71 


0.41 


4.3 


NGC 5055 


0.89 


8.9 


1.00 


5.1 


1.58 


1.00 


1.41 


4.06 


2.0 


NGC 7331 


0.22 


4.5 


0.50 


6.8 


1.58 


0.79 


1.41 


0.67 


2.0 
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